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Abstract 

Recent advances in random matrix theory have spurred the adoption of eigenvalue-based detection 
techniques for cooperative spectrum sensing in cognitive radio. These techniques use the ratio between 
, the largest and the smallest eigenvalues of the received signal covariance matrix to infer the presence 

or absence of the primary signal. The results derived so far are based on asymptotical assumptions, 
due to the difficulties in characterizing the exact eigenvalues ratio distribution. By exploiting a recent 

^ \ 

' result on the limiting distribution of the smallest eigenvalue in complex Wishart matrices, in this paper 

Q>^ , we derive an expression for the limiting eigenvalue ratio distribution, which turns out to be much more 

accurate than the previous approximations also in the non-asymptotical region. This result is then applied 
' to calculate the decision sensing threshold as a function of a target probability of false alarm. Numerical 

_ simulations show that the proposed detection rule provides a substantial improvement compared to the 

other eigenvalue-based algorithms. 

X 

^ • I. Introduction 

Blind detection algorithms, relying on the received signal diversity achieved through multiple antennas, 
user cooperation, or oversampling, have been recently proposed for Cognitive Radio. Most of these 
methods [1], [2] are based on the properties of the eigenvalues of the received signal's covariance matrix 
and use results from random matrix theory (RMT). 

Their main advantage, with respect to classical energy detection (ED) or cyclostationary feature 
detection (CFD) [3], is that they do not require any prior information on the primary signal or on 
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the noise power. Among blind algorithms, the eigenvalue-based approach was shown to outperform ED, 
especially in case of noise uncertainty [2]. 

However, the decision rules of the eigenvalue-based detection schemes proposed so far are based on 
asymptotical approximations, that make them inaccurate in many practical scenarios. Using some recent 
RMT results, in this paper we first derive an analytical expression for the limiting distribution of the ratio 
between the largest and the smallest eigenvalues of the covariance matrix. Then, based on this result, 
we obtain a novel decision rule that outperforms the previously proposed eigenvalue -based detection 
schemes. 

The rest of the paper is organized as follows: Sec. |ll] re views the eigenvalue-based algorithms proposed 
in the literature. Sec. JII] deals with the threshold optimization problem and presents the contribution of 
this paper. Numerical results are presented and discussed in Sec. HV] Sec. |V] contains the conclusions. 

II. Mathematical Background 

A. System Model 

Denote with K the number of collaborating receivers (or antennas) and with the number of samples 
collected by each receiver during the sensing time; let yk{n) be the discrete baseband sample at receiver 
k (k = 1, . . . , K) and time instant n{n = 1, . . . , N). Two hypotheses exist: under (no primary signal: 
the samples contain only noise) yfc(n)|'Ho = vin), where v{n) is circularly symmetric complex Gaussian 
(CSCG) noise with zero mean and variance cr^; under Tii (presence of primary signal) yk{'n)\'Hi = 
hk{n)s{n) + v{n), where s{n) is the primary signal, with E|s(n)p = cr^ / 0, and /ifc(n) is the channel 
between primary source and receiver k at time n. 

Let y(n) = [yi(n) . . . be a K x 1 vector containing K received samples at time n and 

Y = [y(l) . . . y(-/V)] a. K X N matrix containing all the samples received during the sensing period. 
The sample covariance matrix, R(A^) = ;^YY^, converges to R = E[yy^] for N —>■ oo: from the 
eigenvalues of R(A^) it is possible to infer the presence or absence of primary signal. 

B. Previous Results 

Let Xmax and Xmin be the largest and the smallest eigenvalues of R(A^), and Imax and Imin those of the 
normaUzed covariance matrix, defined as R'(A^) = ^R(A^). Under Hq, R'(A^) turns out to be a complex 
white Wishart matrix and, by the Marchenko-Pastur law, the eigenvalue support is finite [4]. Under Hi, the 
covariance matrix belongs to the class of 'spiked population models' and its largest eigenvalue increases 
outside the Marchenko-Pastur support [5]. This property suggests to use T = Imax /Imin = Xmax /Xmin 
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as test statistic for signal detection. Denoting as 7 the decision threshold, the detector decides for if 
r < 7, for Til otherwise. Two approaches to set 7 are proposed in the literature. 

1) Asymptotic Approach [1]: Thanks to the asymptotical properties of Wishart matrices [4], the 
smallest and the largest eigenvalues of R'(A^) under Ho) converge almost surely to 

Imin ^ a = [n'/^ - K^/^y (1) 

Imax ^ b = [n'/^ + K'/^y (2) 

in the limit 

N,K -^00 with K/N c (3) 

where c G (0, 1) is a constant. Under Hi, according to the theory of spiked models, the largest eigenvalue 
converges almost surely to a value b' > b [5]. Based on these results, an asymptotic detection rule was 
proposed in [1] with decision threshold 

b 

las = - (4) 

a 

2) Semi-asymptotic Approach [2]: This approach is based on the use of the recently-found limiting 
distribution of Imax instead of its asymptotical value ([2]). Results from [6] state that under the same 
assumptions (O the random variable 

T Imax b 

J-'max — \jj 
V 

with 



converges in distribution to the Tracy-Widom lawo of order 2. The authors of [2] exploit this result to 
link the decision threshold to the probability of false alarm, defined as 

Pfa = P{T > 7|Ho) (7) 

by using the asymptotical limit ([T]) for the smallest eigenvalue and the Tracy-Widom cumulative distri- 
bution function (CDF) for the largest one. The threshold can be written as: 

Isa = las ■ (1 + ^^^^^0W^PtU^ - Pf^)) 

where F^^2{y) the inverse Tracy-Widom CDF of order 2. 

'The Tracy-Widom distribution was defined in [7] as: 
Ftw2{s) = exp ^— ^ s)q'^ {x)dx^, wliere q{s) is the solution of the Painleve II differential equation q"{s) = sq{s) + 

2q'^{s) satisfying the condition q{s) ~ —Ai{s) (the Airy function) for s +00. For its importance in RMT this distribution 
has been extensively studied and tabulated; a Matlab routine to compute is available at [8]. 
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III. Eigenvalue Ratio Distribution and New Detection Threshold 

The asymptotic approach (Sec. III-Bll ) uses Umiting approximations, valid for very large N and K. 
In practical conditions, that may be characterized by small number of observations due to time-varying 
channel and/or detection in the shortest possible time, the asymptotic threshold turns out to be very 
unbalanced with respect to the actual eigenvalue ratio distribution (see next section. Fig. [fl. In addition, 
this approach does not allow to tune the threshold as a function of a target Pfa- The semi-asymptotic 
approach (Sec. III-B21 ) allows such a control, but it is still based on the asymptotical limit for the smallest 
eigenvalue and it becomes inaccurate when decreases. 

Recently, Feldheim and Sodin [9] found that the smallest eigenvalue also converges to to the Tracy- 
Widom distribution as K,N —>^ oo, up to a proper rescaling factor. Thus, the random variable: 

T . — ~ " fQ\ 

^J' 

converges in distribution to the Tracy-Widom law of order 2, with: 

/i = - (^-1/2 _ ^-1/2^ V3 ^^^^ 

As a consequence of (O, fj, is always negative in the considered range of c. Now, the test statistic T may 
be written as: 

rp ^max ^L'max ~\~ b (11) 

Denote with fi^^^{z) and fi^^.^{z), respectively, the limiting probability density functions (PDFs) of 
the numerator and the denominator of T for K, N ^ oo. From (|5]l and these PDFs may be expressed 
through a linear random variable transformation of the second-order Tracy-Widom PDF, fTW2 (x) : 

lln^aA^) = -fTW2(- ^ 



V 1^ 



(12) 



and, recalling that jj, < 0: 



Jl^J^) = J^fTW2 (^) = --fTW2 ('-^) (13) 

Finally, assuming fi^^^^^ {z) and f^^^^^ {z) as independent (which is reasonable for limiting distributions, 
with the size of R'{N) tending to infinity) and applying the ratio distribution formula [10], we can write 
the PDF of T as: 



oo 
■foo 



\x\fi,,i^{tx,x)dx 
xfi,{tx)fi^^{x)dx 



I{t>l} 



■ I{t>i} (14) 
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Eigenvalue ratio CDF - N=1 000 K=50 
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Fig. 1. Eigenvalue ratio CDF obtained using the novel ratio-based approach vs. empirical CDF, asymptotical approach, and 
semi-asymptotic approach. A'^ = 1000, K — 50. 



where the lower integration limit has been changed to instead of — oo, since the co variance matrix is 
positive-semidefinite therefore all the eigenvalues are non-negative; /{.} is an indicator function, with the 
condition t > 1 to preserve the order of the eigenvalues (h > Ik)- 

Given this new result, we can now introduce a sensing algorithm based on this limiting eigenvalue 
ratio distribution. Let be the cumulative density function (CDF) corresponding to (fT4l ). From ([7]), 

the false alarm probabiUty is Pfa = 1 — Fxi'y), for large N and K; hence, we derive the the novel 
decision threshold as a function of the false-alarm probability: 

Jrd = FT\l-Pfa) (15) 

In practical applications the values of F^^{.), evaluated numerically off-line, can be stored in a look-up 
table and then used by the receiver to set the proper threshold as a function of N, K, and the target Pfa- 
(Note that a look-up table or a similar approach is also needed for implementing since Frp^2 does 
not have a closed-form expression). 

IV. Numerical Results 

Fig. [T] represents the eigenvalue ratio CDF resulting from the novel analytical approach and compares 
it to the empirical distribution, computed by Monte-Carlo simulation, and to those obtained from the two 
approaches of Sec. III-BI The number of samples was set to = 1000 and the number of cooperating 
receivers to = 50. The novel analytical CDF matches with the empirical data, whereas the asymptotic 
one (which is simply a step function) and the semi-asymptotic one are very unbalanced because the 
considered parameters (N = 1000 samples and K = 50 receivers), although large, are still far from the 
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Completnentary-ROC: N=1000 K=50 SNR=-21clB 
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Fig. 2. Complementary ROC: the novel ratio-based approach vs. asymptotical approach, semi-asymptotic approach, and energy 
detection. = 1000, K = 50, SNR = -21dB. 

asymptotical region. From the detector's point of view, this means that neither the asymptotic nor the 
semi-asymptotic approach allow to set the decision threshold correctly according to the target Pja- 

Fig. |2] provides a performance comparison of the considered eigenvalue-based detectors, plus the 
traditional energy detector using a cooperative equal gain combining scheme [11]. This type of graph, 
commonly used for signal detection and called Complementary-ROC (Receiver Operating Characteristics), 
represents the achievable probability of missed detection Pmd = P{T < tI^i) vs. the target Pja- The 
simulation parameters are again N = 1000 and K = 50; the average signal-to-noise ratio under Hi, 
defined as SNR = -li^^ with ||hf = Ef=il^fc|^ is equal to -21 dB. Such low values of SNR 
are typically used to evaluate detectors in critical conditions (e.g., in the case of "hidden node"). For 
energy detection, a noise uncertainty of 0.25 dB is assumed, whereas the eigenvalue-based algorithms are 
insensitive to the noise power uncertainty. The ROC plot shows that the novel ratio-distribution threshold 
provides lower probabilities of missed detection than the other approaches for any given probability of 
false alarm. Since the new algorithm uses a nearly-exact distribution, it allows to choose the lowest 
possible threshold for a given target Pfa, i-c, to obtain the minimum value of Pmd- 

For instance, given a target Pfa of 10^^, the novel approach provides a Pmd of 1.0 • 10^^, while the 
semi-asymptotic approach would give 6.5 • 10^^. The asymptotical approach, as previously mentioned, 
does not allow any control of Pmd vs. Pfa since the threshold is fixed. The pair of {Pfa, Pmd) it achieves 
is represented by a dot in the figure, at (4 • 10~^, 1.15 • 10^^); this value of Pmd = 1-15 • 10~^ is a lower 
bound that cannot be improved regardless of the target Pfa, as highlighted by the straight dashed line. 
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V. Conclusion 

In this paper an expression for the Umiting eigenvalue ratio distribution in Wishart matrices has been 
derived and it has been appUed to the problem of signal detection in cognitive radio. The analytical 
distribution has been shown to be consistent with the empirical data and, for this reason, the novel 
detection rule clearly outperforms the previously proposed ones especially for realistic numbers of sensing 
samples and cooperative receivers. 
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